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ABSTRACT 

This article gives an analysis of the behavior of polarizing grids and reflecting polar- 
izers by solving Maxwell’s equations, for arbitrary angles of incidence and grid rotation, 
for cases where the excitation is provided by an incident plane wave or a beam of ra- 
diation. The scattering and impedance matrix representations are derived and used to 
solve more complicated configurations of grid assemblies. The results are also compared 
with data obtained in the calibration of reflecting polarizers at the Owens Valley Radio 
Observatory (OVRO). From these analysis, we propose a method for choosing the opti- 
mum grid parameters (wire radius and spacing). We also provide a study of the effects 
of two types of errors (in wire separation and radius size) that can be introduced in the 
fabrication of a grid. 

Subject headings: instrumentation: polarimeters — techniques: polarimetric — tele- 
scopes 


1. Introduction. 

The literature on wire grids is abundant and they have been studied with different techniques 
and for numerous applications. Most of the analysis were however restricted to special cases of 
incident field and grid orientations. The more general and arbitrary situation seems to have been 
first studied by Wait (see Wait (1955a) and Larsen (1962)). This problem is addressed again in 
this paper and follows a line of analysis fairly similar to the one used by Wait. Our treatment is, 
however, more general in that we do not assume that the wires of the grid are induced with only 
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a longitudinal current; we will indeed show that an azimuthal component is also present. We also 
solve for the induced current by considering the tangential components of both the electric and 
magnetic fields at the surface of the wires. 

This analysis is carried out in the next two sections and will serve as our basis for the treatment 
of the reflecting polarizer (section 4) and the introduction of the scattering and impedance matrix 
representations for a grid (section 3.3) which will in turn enable us to briefly discuss more compli- 
cated systems. These matrices will be particularly useful in allowing us to define what will be called 
the principal axes of a grid. These are two orthogonal and independent directions of polarization 
in the plane of the incident radiation along which an arbitrary electric field can be decomposed 
and shown to scatter without cross-polarization. With this representation at hand, it will then be 
possible to derive a set of optimal parameters (wire radius and spacing) to be used in the selection of 
a grid. We will also present an analysis of the effects of random errors that can be introduced in the 
fabrication of grids, the results obtained will then be compared to experimental results previously 
published by Shapiro & Bloemhof (1990). 

The last section will be dedicated to the study of the more subtle impacts that the nature of 
the incoming radiation can have on the response of a grid assembly such as a reflecting polarizer 
(section 4). Although limited to this particular case, our discussion could possibly apply to other 
types of instruments. We have also included at the end (Appendix II) a list of the symbols used in 
the different equations. 


2. The case of a single wire. 

Before trying to solve the problem of the grid or the reflecting polarizer, it is preferable to 
study the case of a single conducting wire. It will serve as the basis for our studies of the more 
complicated cases to follow in subsequent sections. 

Let's suppose that a wire of radius a is oriented, as depicted in Figure 1, parallel to the x-axis 
at y — y Q , z = and that it is subjected to an incident plane wave E*(r) of arbitrary direction and 
polarization: 


Ej(r) = E„(a'e x + 0'e y + 7 'e 2 )exp(-j(k ■ r - wt)) (1) 

with 


k = k(ae x + (3e y 4- 7 e 2 ) 

and where, of course, the following conditions of normalization and orthogonality apply: a 2 + /3 2 + 
-y 2 = a ' 2 4- ; 7' 2 4- ~j rI = 1 and aa' + fif}' + 77' = 0. Using the coordinate system depicted in Figure 1 
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we have a = sin(xi) sin({p 5 ), /3 = sin^) cos(tp g ) and 7 = cos(xi) where \i ls the angle of incidence 
and cp g the angle of grid rotation. 

In everything that follows, we will drop the exp(jn?£) term and assume it to be implicit in the 
equations. We will also suppose that the wire is of infinite length and made of a good conducting 
material of conductivity a such that any current flowing through it can be accurately represented 
by a surface current vector K. This quantity is related to the current density J(r) as follows: 

J(r) = K S(p - a) cxp(-jk • r) (2) 

where 

K = K x e x + K e e e 

and y — y Q = pcos(8) : z — z Q = psin(0). 

Before we solve for the scattered fields, it is to our advantage to note that for the case considered 
here (i.e., thin wire with an approximate solution involving no angular mode dependency), the 
problem can be broken in two parts or modes. The mode where the electrical field is parallel to 
the plane defined by and k (the transverse magnetic or TM-mode) is related to the presence of 
K x while another, where the magnetic field is parallel to this same plane, the transverse electric or 
TE-mode is related to K 6 . The analysis will, therefore, be facilitated with the use of the two vector 
potentials A s and for the scattered fields (Balanis 1989). 

The TM-mode can be analyzed using the vector potential A s , in the Lorentz gauge, with 
F s = 0. The needed equations are: 


-Mr) 

= /j(r') eXP( “ ,tfi) <(V 

47 1 j R 

(3) 

E s (r) 

= ^V(V-A s (r))-jwA,(r) 

= — VxA s (r) 

(4) 

H s (r) 

(5) 


Mo 


with 

R 2 = (x - x') 2 + (y- y’f + (z - z') 2 . 

Since we are concerned here with the longitudinal component of the surface current density, 
we only need to consider the A x component of the vector potential (i.e., we set A p = A 6 — 0). 
Equation (3) can be solved exactly when K is expanded with a Fourier series, but in cases where 
the wavelength of the incident wave is much larger than the wire radius it can be shown that: 
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A x ( r) = ^K x H [ /\k'p) exp(-jV) (6) 

with fc' = fcv'l - <?, V 5 = fc ( a;r + PVo + 7z 0 ) and where H^\x) is Hankel’s function of the second 
kind of order n. 

On the other hand, it is advantageous to study the TE-mode with the vector potential F s , in 
the appropriate gauge, with A, = 0 (Balanis 1989). To do so, we will not consider the effect of the 
current density (more precisely its azimuthal component) but that of the magnetization vector M 
that it induces. The relevant equations are now: 


J( r ) 

P.(r) 

= V x M(r) 

_ Woeo f M(x >) a Vl-j kR ) i p r J 
4tt J *'*' R 

(7) 

(8) 

B.(r) 

= --VxF.(r) 

(9) 

H,(r) 

= — V(V-F s (r))-jwF s (r) 

(10) 


with R as defined above. 

Since we are now concerned with the azimuthal component of the surface current density, we 
only need to consider the M x , F x components of the magnetization and vector potential (i.e., we 
have M p = M 9 = Ff = F* = 0). Again, in cases where the wavelength of the incident wave is 
much larger than the wire radius it can be shown that: 

Ff( r) = 7rWti °^ -— K e H { 0 2) (k'p) cxp(-jy>) . (11) 

It is now straightforward to calculate the scattered fields by combining the solution obtained 
for each mode (using equations (4), (5) and (6) for the TM-mode and equations (9), (10) and (11) 
for the TE-mode): 


E p (r) - -jay/ I - ot 2 F K x H { /(k'p) exp(-j<p) (12) 

E 0 S ( r) = -y/l^cPF y K e H[ 2 \k'p)exp(-j<p) (13) 

E x ( r) = - (l - a 2 ) FK x H^\k'p)exp(-j(p) (14) 

Jff(r) = a F^ K 6 H ( /\k'p) cxp(-jy) (15) 

Z 0 l 
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Hf{ r) = -j Vl „ Q2 F K x (k‘ p) exp(-jip) (16) 

Hf(r) = -j ^ ~“ 2) F y K e H$\k'p) cxp(-jV) (17) 

where F = and Z 0 = is the impedance of free space. Note that although equations 

( 12)- (17) represent the scattered field, the components of surface current density that are included 
in these equations axe that of the, yet undetermined, total surface current density which we are now 
in a position to evaluate. 

In order to do so, we must first express the incident plane wave in the appropriate coordinate 
system. This can be done by first using the following expression: 


oo 

e*p(~jk(Py + 7*)) = E (-j) nj n(k'p)e\p(jnO') (18) 

fl — — OO 


with 6 f = 9 — arctan ^ j , J n (x) the Bessel function of order n and by again splitting the incident 
field in the two modes defined earlier (van de Hulst 1957; Balanis 1989). This enables us to express 
the plane wave in cylindrical coordinates and match the fields with the usual boundary conditions 
for their tangential components at the surface of the wire. For the TM-mode the condition is: 

Ef + Ef = Z s (Hf + Hf) . (19) 

Within the order of precision used for our analysis (A a) and considering a solution with no 
angular dependency, it can be shown that: 


Ef ~ a'E 0 exp (-jk (ax + 7 z a )) 

Hi ^ ja' y ■ y ex P ( ~jk (ax + 7 z a ) ) . 

For the TE-mode we have: 

Ef + Ef = - Z s (Hf + H x s ) (20) 

with: 


ka 

Ef ~ -j ( 7 'P - & 7 ) E 0 — exp (-jk (ax + 7 z a )) 
Hf ^ (j'P - fi'l) rr exp(-jk (ax + 7 z a )) . 

^ O 
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In equations (19) and (20), Z s = (1 4- j)y * s ^e sur f ace impedance of the wire (Jackson 
1962). It is to be noted that for wires of small radius, relative to the wavelength, the boundary 
conditions (19) and (20) along with the equation for Z s represent approximations that are only valid 
in the lowest mode and for a sufficiently good conductor. A more rigorous treatment shows that 
these equations will be modified in the more general case (Wait 1979; Bouche, Molinet & Mittra 
1997). But for the purpose of our analysis, the approximation used here is adequate. 

When solving these two sets of equations we find the following expressions for the components 
of the total surface current densities: 


E 0 

“'(I ->£•*) 

(21) 

F 

(1 - a 2 ) H ( a\k’a) 1 - a 2 Hf\k'a) 

E 0 


(22) 

F ' 

VI - a 2 H?{k'a) + (1 - a 2 ) H^(k'p) 


These last two equations can be inserted in equations ( 12)-( 1 7) to calculate the value of the 
fields at any point exterior to the wire. For a good conductor the internal fields are practically 
nonexistent. Equations (21) and (22) are in agreement with the results presented in Balanis (1989, 
see ch. 11) for the case of normal incidence and a perfectly conducting wire. 


3. The polarizing grid. 

3.1. Analysis. 

With the solution for a single wire in hand, the problem of a configuration of an infinite number 
of wires of infinite length separated by a distance d is simplified if one realizes that every wire will 
be induced with the same surface current K. The only difference will be a phase term in the current 
density J(r), given by equation (2), which depends on the position of the wire along the y-axis. The 
same thing can be said for the scattered fields from any given wire, one only has to replace y Q by 
nd in equations (12)-(17), where n is an integer that determines the position of the wire. 

If the scattered fields are now just the sum of all the different scattered fields from the individual 
wires, care must however be taken in evaluating the surface current. First, when one matches the 
boundary conditions it must be done simultaneously at the surface of every wire. However, since 
we are dealing with an infinite number of infinitely long wires subjected to the same incident plane 
wave, it turns out that it is sufficient to do so for only one of the wires. If the boundary conditions 
are matched for one wire they will be for all. We have chosen for our calculations the “center” wire 
at n = 0. Second, to match the boundary' conditions we must express the scattered fields of each 
and every wire in a cylindrical coordinate system centered on the position of this “center” wire. 



When this is done, we find the following expressions for the components of the induced total 
surface current density: 


with 


and 


K x 

K° 


E 0 f N x 
F a A x 

• — ° ( f ft ft 9 \ 


N x 

A* 

Ng 

A<? 


= 1 -3 


. Z s ka 


Z 0 2 


= (1 - a 2 ) Si - jf 1 • V 1 - a?H™(k'a) 


i , ■ 2 

3 Z 0 ' ka 


Vl -a 2 i/{ 2) (A:'a)+i|^ ■ (l - a 2 ) 5! 


S\ = + 2 ^2 H™(k'nd) cos (kfind) . 


n= 1 


(23) 

(24) 


(25) 

(26) 

(27) 

(28) 


(29) 


We will give in section 3.4 adequate approximations for A x and A# that will greatly simplify 
the evaluation of the reflection and transmission coefficients which are soon to follow. 

By using the appropriate expansions for series of Hankel’s functions we can write down the 
expressions for the components of the total electric field far away from the grid: 


m r) 
E y T (r) 


(l — X F 

= Q'£^ 0 exp(-jk • r) -K x expi-jk^z - z 0 \) exp(-jtp) 

7 7 Tu 


= p'E 0 exp(-jk • r) 


+- 


\F 


E r( r ) = l'E 0 exp(-jk ■ r) 


K e z ~ z q 


Z - z 0 


+ 


A F 

nd 


aK 1 


z - .0 ka g 

\z-z 0 \ 3 7 2 


exp{-jk^\z - z 0 \) exp(-jip) 
exp(-jky\z - z 0 \)exp(-j<p) 


where p = k(ax + fly + 7 z 0 ). From these it is now straightforward to get the reflection and 
transmission coefficients (normalized to E 0 ) in the far-field: 
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R x 

R y 


R z 

rpX 

J'V 


T z 



A 

ird 

(! « 2 ) R x 

7 


A 

7T d 

* 

H 

Hr 

>! 

<*> 

1 1 


A 

■nd 

L K r + /_ ^ K , 
7 z J 


at' + R x 


&>Jr Y 0 Vd 

, F A 

V H — 

1 E 0 7 rd 


Ck 


(30) 

(31) 

(32) 

(33) 

(34) 

(35) 


where we have set z 0 = 0 for simplicity. 

Equations (30)-(35) along with (23)-(24) arc the solution to the polarizing grid problem for 
cases where it is assumed that k'a <C 1 and a<rf. 

For predictions of measurements made in the laboratory, one merely has to transform these 
coefficients to the laboratory coordinate system. If we adopt for this system the coordinates of the 
incident/transmitted («, v,w) and reflected (u',v',w') plane waves defined in Figures 1 and 2, the 
last system of equations is simplified to: 


R u 

R v ’ 

R w ' 

rjvii 

rjiV 

/jfW 


F 

E~o 

F 

~Eo 


7T d 

X_ 

7T d 


1 

7 y/l — 7 2 
1 

7\/l — 7 2 


l3K x -j ai ~K e 
07 K x + 3 P^K e 


1 


a 




E 0 nd 7 -y/l — 7 2 

A 1 


Eo 7 T d 7^/1 - 7 2 


pK x +jaj^-K e 
al K x -jp^K 6 


with a" and /?" related to the incident field by : 

Ej(r) = E 0 (o"e u + 0"e v ) exp(-jkw) . 


(36) 

(37) 

(38) 

(39) 

(40) 

(41) 


As can be seen, the reflected and transmitted fields have no component along their respective 
direction of propagation as is required for the propagation of plane waves in free space. 
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3.2. Effects of grid imperfections. 

So far we have assumed that there were no imperfections in the construction of the grid, 
obviously (and unfortunately) such is not the case in a realistic situation. It would be instructive 
if we could calculate the effects of errors that are likely to be introduced in the fabrication process. 
In this section we will provide expressions that will allow us to evaluate changes in the reflection 
and transmission coefficients induced by two possible imperfections: random errors in wire spacing 
and random variations in the size of the wire radius. 


3.2.1. Random errors in wire spacing. 

It is our experience that some of the commercially available grids when observed under a 
microscope show some defects in their assembly. Visually, the most obvious manifestation of this is 
inconsistency in the spacing between wires. In order to calculate the effect of these errors we have 
to go back to the discussion of section 3.1 that guided us into the evaluation of the induced current 
on the wires. Since we can no longer assume that the wires are evenly spaced, we must now realize 
that they will in general have different values for the current and fields on their surface. This wall 
be made more apparent if we write down the expression for the ^-component of the electric field on 
the surface of the “center” wire: 


oo 

E g( a ) ~ ~ (! “ q2 ) f exp{~jk(ax + 720 )) K% (£)<?„ (0 ( 42 ) 

n=- 00 

with: 


Gn (0 = 


Hl 2) (k'a)exp(-jkp£ 0 ) ,n = 0 

Hl 2) (k 1 |nd + £ n |) exp(—jkP (nd + £„)) ,n ± 0. 


(43) 


K* is the induced surface current on wire n and the £ n axe statistically independent random 
errors in the positioning of the wires. Now, if E {x} stands for the expected value of x and if we 
suppose that the errors have a zero mean, we can write: 


EiCm} 

= 0 

, r = 1, 3,5 . . . 

(44) 

E{C m } 

= e {&} = £{?} 

,Vm,n 

(45) 

E{CmC} 

- £{&}£{£} 

, m ^ n 

(46) 

i m K x 


, Vm,n . 

(47) 


The first equation is deduced from the supposed evenness of the probability density function of 
the errors, the second states that their statistics are the same across the grid and the third expresses 
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their statistical independence. The last of these equations arises from the fact that if we were to 
test a large number of similar grids, every wire would exhibit the same average value £ m K x for any 
induced surface current moment (independent of its position n). 

We will not go into the details of the calculations as they are somewhat lengthy, but it can be 
shown that if we apply this last set of equations and expand (£) and G n (£) with their Taylor 
series around = 0 while solving for the boundary conditions, we can find an expression (valid to 
the second order in f) for the average longitudinal surface current: 


K x ~ K x 



(1 - a 2 ) 


\dGm] 

to 

1 


' £m — 0 j 


where K x is the current density induced on the wires of a perfect grid and is given by equation (23). 
One sees that the errors bring a perturbation which is proportional to their common variance. 

If the same approach is used to calculate the effect of such random errors on the value of the 
azimuthal surface current density K one finds that it remains unaffected: 


K e ~ K 6 


with K e given by equation (24). 

From this we could then proceed and calculate the expected value of the reflection and trans- 
mission coefficients by evaluating equation (42) (and the corresponding equations for the y and 2 
directions) in the far field, when this is accomplished we find that the coefficients have exactly the 
same form as shown in equations (30)-(35) (or (36)-(41)). We then merely have to replace K x and 
K e by K x and K 6 respectively. 


3.2.2. Random errors in the wire radius (wire to wire). 

Another type of error which can be analyzed is one concerning the random variation in the 
size of the wire radius, which we will denote by the letter r/. More explicitly, we are considering 
differences between wires and not variations along a single wire; we assume the diameter of a wire 
to be constant but somewhat uncertain in its value. This is the kind phenomenon that could occur 
if the wires were stretched with slightly different tensions when installed or perhaps also in cases 
where the wires have a finite ellipticity and are rotated between rows. We can proceed in the same 
manner as we did in the last section for the analysis of the boundary conditions and the fields away 
from the grid. When this is done we get: 
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K* ~ K x 1 
K° ~ K 6 1 


E{V 2 } 


£“=- ocG n ( 0 ) 


f 1 d 2 G 0 

(i - or 2 ) 

' OGo ■ 

n 

{ 2 drj 0 2 

Ax 

. dT )0. 

' Vo=0_ 


£{t/ 2 } [1c> 2 Qo ldQp 

Qo (0) [ 2 <fy 0 2 O dr]o A# [ ctyo J 


9Qo 


} 

' V 0=0 


where r) n is the random error in the size of the radius of wire n, G n is given by equation (43) 
(with £ = 0 and a replaced by a 4* r j), A x and A# by (26) and (28) respectively and Qo(r)o) — 
h[^ ( k f (ao + tjq)). Again the expected value of the different coefficients can be obtained by replacing 
the current components K x and K e by and K 6 in equations (30)-(35) (or (36)-(41)). It will 
also be noted that the errors contain a perturbation term which is proportional to their common 
variance. 


3.2.3. Predictions and comparison with experiments. 

Now that we have derived the equations for the reflection and transmission coefficients it would 
be interesting to compare the predictions that our model makes with experimental data. Although 
we have independently treated the two types of errors, it is nevertheless obvious that within the 
limit of precision of our analysis (small errors) that they can both be simultaneously added in 
the expressions for the reflection and transmission coefficients. Doing so would in principle allow 
us to compare theory and experiments as actual grids are liable to exhibit both kinds of defects. 
This also suggests though that it might be impossible to separate the effects of both errors in 
measurements. It turns out, however, that the perturbations caused by the errors in the size of the 
wire are predicted by our model to be smaller than those caused by the errors of the other type 
(for equivalent error amplitudes), and we neglect them in the following comparison of theory and 
measured grid properties. 

Shapiro & Bloemhof (1990) have published measurements of the unwanted cross-polarized 
transmittance through three grids on which they had purposely introduced random errors in the 
wire positioning. They quoted the errors in term of the random variation in the distance between 
wires (pitch) with amplitudes of 7%, 23% and 52% of the mean wire separation (aimed at 108 /im 
with a wire radius of 12.5 /im). We must divide these values by a factor of y/2 in order to relate 
them to our errors £ n since we have defined these as pertaining to the absolute position of the 
wires. Figure 3 shows a comparison of our model’s predictions with their measurements done for 
cases where the incoming field is at normal incidence to the grid and polarized parallel to the 
wire orientation. Although the agreement is not perfect, the outcome is very satisfactory as the 
theoretical curves exhibit the right behavior with frequency and error amplitude. 
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3.3. The scattering matrix and the impedance model. 

3.3.1. The scattering matrix and the principal axes of a grid. 

The relationship between the reflection and transmission coefficients in equations (36)-(41) is 
reminiscent of what is often encountered in microwave engineering in the analysis of systems that 
can be accurately dealt with using a lumped-elements model. With this in mind, it is tempting to 
consider any problem involving a polarizing grid by treating the different components as lumped 
and interconnected through a transmission line of characteristic impedance Z 0 (Lamb 1997). We 
can then go ahead and model the grid as a 4-port device since the reflection and transmission 
coefficients given by the aforementioned set of equations provides us with the scattering parameters 
at each port. 

In this context, it is more convenient to work with a single coordinate system (u,v,w) (see 
Figures 1 and 2) for both the incident/transmitted and reflected plane waves since we can assume 
that their propagation is done along the same transmission line (it is however understood that, in 
reality, away from normal incidence the transmitted and reflected waves travel along different axes). 
We therefore assume that the incident/transmitted fields travel along the w - axis (with the u-axis 
vertical and the v-axis horizontal) and the reflected fields along the negative w-axis as seen from a 
given side of the grid. 

Since there are two possible independent states of polarization (with the field aligned along the 
u or u-axes), where the waves can travel either toward or away from the grid, we need two ports 
on each side of the grid. So for example, if the incident wave on a given port has an electric field 
polarized along a given axis we can define 4 scattering parameters: one for the reflected signal at 
the input port and three for the transmissions to the other ports. The same thing can be done for 
every port leading to a total of 16 scattering parameters. 

In what follows, a scattering parameter s mn is defined with the 3 ports m n terminated with 
the line characteristic impedance Z 0 . Also, each port n has two signals: an incoming signal £„ and 
an outgoing signal E~;n = 1, 2 (3, 4 on the other side of the grid) refer to polarization along the 
u and v axes respectively. The scattering matrix relates the different signals as follows: 


'£r' 


sn 

«12 

3 13 

Si4 

E; 


521 

S 22 

323 

324 

Ez 


531 

3 32 

333 

334 

. E 4 . 


_ 541 

342 

343 

344 


Et 

E} 

E t 

E* j 


(48) 


The elements of the matrix can be directly evaluated from equations (36)-(41) and shown to 


be: 
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with: 


R uu 

1 R uv 

rpuu 

rpuV 

R uv 

R vv 

rpuV 

rpVV 

'J'UU 

rpuV 

R uu 

R uv 

rpUV 

'J'VV 

R uv 

R vv 


R uu 

R vv 

R uv 

rpUU 

rpVV 

rj^UV 


7 r d 7 (1 - 7 2 ^ 

) 


A 1 


a 2 

7r d 7 ( 1 — 7 2 ^ 

) 

A a/3 

\N X 

7 rd (1 - 7 2 ) 

Ax 

A 1 



7 rd 7 ( 1 — 7 2 

) [ 

A 1 


. r 

nd 7 ( 1 - 7 2 ) [ 

A a(3 

' N x 




2 A* 


+ 


2 A e \ 
a N * 


2 2 ho> N$ 

K a 7 — • — 

A x 2 A#_ 


9 9 -^Vr r.9 rvi* 

' 2 ~ 2 — + 0 1 — 


a 7 


/I U/J 

7T</ (1 — 7 2 ) 


A x 
feo TV,} 
T' Ag. 


2 

,2 JVfl ' 

T'Ao. 


where N x , A x , Ng and A# are given by equations (25), (26), (27) and (28) respectively. 


( 49 ) 


(50) 

(51) 

(52) 

(53) 

(54) 

(55) 


We can go one step further and render things considerably simpler if we make a change of 
coordinates and use the following as eigenvectors instead of e u and e v : 


_ /3e u + aryey 
Pl “ y//'i 2 + a 2 7 2 

C*7^u “1“ ft&v 

^ y/ fi 2 + a 2 7 2 

From now on wc will refer to these as the principal axes of the grid (for reasons that will soon 
become apparent). A close examination of the first of these two equations shows that pi is parallel 
to the projection of the direction of the wires in the plane of the incident field. The matrix S then 
takes a simpler form (we also interchange the second row with the third and the second column 
with the third): 


S - 


Ry Tj| 0 0 

7]| i?|| 0 0 

0 0 R± T ± 

0 0 r x il| 


(56) 
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where 


R = A (l-« 2 ) Nx 
^ 7 rd 7 A x 

= (1-a 2 ) a Np 

X 7 d A# 

T,| = l+i?n 

T x = l-i?x- 


(57) 

(58) 

(59) 

(60) 


We then have a further simplification in the modeling of the grid, evidently equations (57)-(60) 
represent the reflection and transmission coefficients along the two principal axes. 

This last representation has the advantage of simplifying calculations since it allows us to de- 
compose any incident field into two non-interacting components, one along each one of the principal 
axes. That is, a field polarized along one of the principal axes scatters only in this same polarization 
state (as can be deduced from the block-diagonal form of equation (56)). It is also interesting to 
note that even though we have defined the principal axes within the framework of our approxima- 
tion of the grid (fc'o « 1 and a < d), the result obtained here still holds in the general case (see 
Appendix I for a proof). This implies that for the case where one wishes to use a different approach 
to solve (numerically or otherwise) the scattering off a grid of arbitrary characteristics, it will always 
be possible to split the incoming field along the principal axes therefore avoiding cross-polarization 
terms and greatly simplifying the solution. 


3.3.2. The impedance model. 

It seems reasonable to think that a grid could also be modeled with another representation 
where the scattering matrix is replaced by an impedance matrix which contains the same number of 
elements since as before, the grid is still treated as a 4-port device. In this case however, the matrix 
relates the total voltages (electric fields) and currents (magnetic fields) between each and every 
port (Collin 1992). The scattering matrix formulation follows more naturally from our analysis and 
has the advantage of dealing with quantities (reflection and transmission coefficients) which are 
directly measurable whereas impedances are not (at least at the wavelengths considered here). The 
impedance model has however received a great deal of attention in the literature and often seems 
to be the way in which polarizing grids are characterized (Wait 1954, 1955a; Larsen 1962). 

Taking advantage of the principal axes representation, it is possible to treat each two-dimensional 
block of the scattering matrix (equation (56)) separately. It can be shown that the impedance matrix 
Z* corresponding to a given block Sj, can be expressed as: 


Z 6 = Z 0 -(I±S ( ,)-(I=fS 6 )- 1 


(61) 
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where I is the unit matrix and the upper and lower signs correspond respectively to the upper 
left and lower right blocks of the scattering matrix (equation (56)). Applying this last equation to 
equation (56) we get: 


Z p Z p 0 0 

0 0 z n -z n 

.0 0 -Zn Z n 

with: 

z = _Zo 1 + fln 

p 2 ' *11 
7 Zo 1 — R± 

" ~ ~2 HZ ‘ 

It follows quite naturally from equations (57)-(60) that we could have defined two impedances 
Z || and Z ± : 


Z i 


Z 0 ■ 


1 + *ll 
l"*!i 


Z„- 


1 - Rx 

1 4" R± 


One can easily verify that and Z± are respectively equal to Z p and Z n placed in parallel 
to the characteristic impedance Z Q . We therefore see that the impedance matrix gives the “actual” 
impedance of the grid along each of the principal axes whereas the scattering matrix includes, as 
should be expected, the contribution of the loads of characteristic impedance Z Q which is assumed 
to be connected to the appropriate ports when defining its parameters. 


3.4. Approximations and selection of a grid. 

We will now study more closely our simpler equations (57) and (58) for the reflection coefficients 
i2|| and i?j_and try to find relations that will allow us to find a set of optimum parameters for the 
selection of a grid. But before we do so, it will be to our advantage to approximate the expressions 
for A x and (equations (26) and (28)). 

So if we limit ourselves to situations where d A, a A and Z s <C Z Q (good conducting wires) 
and use the proper expansion for Hankel’s functions and series of Hankel’s functions applicable in 
such cases (small arguments) we find: 
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where ~ 0.577215 is Euler’s constant. In equation (63) we have kept things to the lowest order 
possible and we did not expand Z s ; the same is not true for equation (62) for reasons that we shall 
encounter shortly. 

We turn now to the problem of selecting the right parameters for a grid. If we decompose 

a given incident field into two components along the principal axes pi and p 2 (see section 3.3.1), 

a perfect grid would completely reflect the first of these and transmit the second (i?j| = -1 and 

2 

Tj_ — 1). As we will soon see, the coefficient of reflection R± is proportional to when Z s 0 
and is therefore a very small quantity for the cases considered here and we will not worry about it 
anymore (i.e., Tj_ is nearly equal to unity). The condition of total reflection will dictate our choice 
for the parameters of the grid. A close study of equation (57) tells us that in order to achieve perfect 
reflection we must simultaneously satisfy the following relations for the real and imaginary parts of 
A x (for what follows we assume N x ~ 1, see equation (25)): 


Re {A*} 
Im {A x } 


A (l-a 2 ) 

7T d 7 

0 . 


Solving for these we then get: 


a 

d 


L(l-a 2 )V<rZj 

27TO . 


(64) 

(65) 


Had we kept equation (62) to the lowest order, we would have been unable to specify an 
optimum value for the wire radius but only the relation that binds d to a. It is also of interest to 
note that for a given wavelength, the finite size of the wire radius is, to this level of approximation, 
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dictated by the conductivity o\ if we let o — ► oc then there is no restriction on the smallness of the 
radius. 

In a quantitative example to demonstrate the values that can be expected for a and d y assume 
that we are working at normal incidence at a wavelength of 1 mm with a grid made of copper 
(a = 5.8 x 10 7 f2 _1 m _1 ). Using these, we obtain a ~ 11 /im and d ~ 70 /im. 

At this point it is appropriate to discuss the implications of the two assumptions we made 
at the beginning concerning the wire radius and spacing, namely that k f a <§C 1 and a«(i. It is 
important to make sure that a given choice of grid parameters are well within the boundaries of 
applicability of our model. As a means of determining these boundaries, we simulated the response 
of grids (and assemblies of grids, see section 4) for different combinations of wire radius and spacing 
and made sure that the results obtained were reliable (for example, it is obviously imperative that 
the magnitude of the reflection and transmission coefficients never exceed unity). As it turns out, 
there is a fairly strong restriction linking the size of the wires and the wavelength, but if one makes 
sure that A > 40 a then one seems to be well within safe modeling conditions, a cannot be too small 
either. However, since for a good conductor (again let’s use copper) the skin depth at 1 mm is on 
the order of 0.1 /im, our assumption of the existence of an idealized surface current is more than 
adequate. It seems that the second restriction concerning the spacing of the wires is not as binding 
as the first one. It is clear that d > 2 a for if not the wires would be touching, but it appears that 
everything is fine for d > 4 a. Our proposed optimized values for the grid are therefore justified. 

It is also appropriate to point out that using equations (62) and (63) for A x and A q (with 
or without the optimal values for a and d given by equations (64) and (65)) along with equations 
(57)- (60) for the reflection and transmission coefficients along the principal axes renders the task of 
calculating the response of a grid a rather simple one. It becomes unnecessary to confront the more 
intimidating representations derived earlier in section 3.1 (compare with equations (23)-(35)). For 
example, to the lowest order, we get for the reflection coefficients: 


i +;¥>»(£) 

. (l - a 2 ) 7 r 2 a 2 

n . 

J A d 


( 66 ) 

(67) 


which are in agreement with known results (Larsen 1962) (more precisely, for the case of normal 
incidence discussed in Larsen (1962), equation (66) reduces to the result presented there whereas 
equation (67) differs by a factor of two or three depending on which approximation it is compared 
to). 
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4. The reflecting polarizer. 

4.1. Analysis. 

From the solution of the polarizing grid it is a somewhat natural extension to consider the 
more complicated problem of the reflecting polarizer. A reflecting polarizer consists of an assembly 
where a polarizing grid, like the one studied in the last section, is followed by a mirror paralleling 
it at some distance z a behind (effectively placing the mirror at z = 0). 

It is appropriate in this case to use the method of images to solve this problem (Wait 1954). 
We then assume that images of both the incident and scattered fields arc emanating from the other 
side of the mirror. This is equivalent to saying that the image world is made of a grid positioned 
at z = -z 0 with an image incident field impinging on it. Assuming that the mirror is made of a 
material of good conductivity, one can write for the image incident field E'(r): 


E'(r) = E 0 {a' m e x + 0, ' m e y - p m e z ) exp(-jk(ax + 0y - -yz )) 

with: 

a' m = tj [a' (a 2 Rtm + P 2 Rte) + ( Rtm ~ Rte )] ( 68 ) 

P' m — ^ _ —r 2 ) [ft (a 2 Rte + ft 2 Rtm) + (Rtm ~ Rte)} (69) 

7 m = i'Rtm ( 70 ) 

where Rte and Rtm are the reflection coefficients of the mirror, with a dependency on the angle of 
incidence, for transverse electric and transverse magnetic modes of incoming radiation respectively 
(Fowles 1975). It is important to note that these tranverse modes of radiation are not the same as 
those introduced in section 2, they are defined here in relation to the plane which is parallel the 
normal vector out of the surface of the mirror (— e z ) and the wave vector k. 

One can go through calculations similar to those carried out in section 3.1 and find the following 
relations between the components of the total surface current densities of the “real” and image grids: 


K' x 

K>0 


7 /^m7 irQ 

YP - Pi 


where K' stands for the surface current of the image grid. 
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From these and by matching the boundary conditions at the grid, it is straightforward, but 
tedious, to solve for the problem. We give here the final results: 


K x 

K e 

R x 

R y 

R z 


Eo 

F 

Eo 

F 


2 ja' sn(/c7h) N x 

(1 — a 2 ) A5i — j^V 1 — a 2 AS 2 
-2j (V/3 ~ £ 7 ) csjkyh) Ne 
Vl -a 2 ES 2 +jf^(l -a 2 )£Si 


exp (-jkyh) 
exp(-jkjh) 


a r 


F 2A (1 - a 2 ) 


E 0 7 id 

of F 2A 
+ “£T ■ J ; 


£*0 7rd [ 7 


j sn(fc 7 /i) exp {jkjh) 
r 

— sn(A:7/t) £ x - ^ cs(A:7h) 


. 2A 
£ 0 ^ 7rd 


/? kfi 

a sn( k-yh) K x H — ■ — cs(k-yh) K e 
7 2 


exp (jk'yh) 
exp (jkjh) 


where 


(71) 

(72) 

(73) 

(74) 

(75) 


sn(a:) = ^ [exp(jz) - r x exp(-jx)] 

cs ( x ) = \ [exp(jx) + r 6 exp(-jx)} 

A5j - H {2 \k'a)-r x H {2 \k'2h) 

OO 

+2 ]T \H (2 \k'nd) - r x H {2) (fcV(nd) 2 + 4/t 2 )] cos (kfind) 

71— 1 

A5 2 = H (2 \k'a)-r x H (2) {k'2h) 

E5i = H(f\k'a) + rgHQ 2 \k'2h) 

OO 

H-2 ^ \H^\k l nd) + {^k f \J (nd) 2 + 4/i 2 j j cos(fc/3nd) 

n= 1 

E5 2 = H[ 2) (k'a) + r e H {2) (k'2h) 


and 

a' 

7mff ~ #»7 
y ^-^'7 ' 

We have also replaced -z G by h (h > 0) so that the distance between the mirror and the grid 
is expressed by a positive quantity. This set of equations along with equations (68)-(70) give us the 
solution of the reflecting polarizer problem for cases where k f a 1 and a d. 
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As was the case for the polarizing grid, if we transform those coefficients to the laboratory 
frame of coordinates ( u ', v ', w') (see Figure 2) we obtain: 


R u ' 

R v ' 

R w ' 


Q r 


3" 

J m 


0 


F 2X j exp ( jk'yh ) 
E 0 7r d 7^/1 - 7 2 
F 1 2A j exp (jhyh) 
E 0 nd 7v /l - 7 2 


0 sn/Zry/i) K x — ay — cs (£7/1) Ff e 
ka 

£*7 sn(fc7/i) Ff 1 + 0 — cs(k'yh) K° 


(76) 

(77) 

(78) 


with a™ and /?((, given by : 


a m = ^ \ a " { a2 RTM + 0 2 Rte) + a/?/?" ( Rtm ~ Rte )] (79) 

/?", = ^ -72 ) ( q2 - r te + 0 2 Rtm) + a0a" {Rtm - Rte)] (80) 

and finally a" and 0" are related to the incident field by: 

E ( (r) = E 0 (a"e u + 0"e v ) e.xp(-jkw) . 


4.2. Solution using the scattering matrix. 

The scattering matrix representation of the polarizing grid gives us the advantage of rendering 
possible the solution of problems that would be otherwise extremely difficult, if not impossible, 
to solve using Maxwell’s equations. For example, a solution of the reflecting polarizer problem is 
straightforward if we “connect” the mirror to ports 3 and 4 of the grid at a distance h behind. Using 
the definitions introduced in the discussion leading to equation (48) we have: 


E 3 = E 3 RTEexp(~jk'y2h) 

E A = E a Rtm exp{-jk-y2h) 

where Rte and Rtm are as defined in the previous section. 

We can then solve for Ey and E^ and find results that are in agreement with those obtained 
in the previous section. 
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4.3. Experimental results. 

Reflecting polarizers like those studied here were tested at the OVRO for polarimetry in the 
wavelength ranges of 1.3 mm as well as successfully used for polarimetry observations at 3 mm 
(Akeson 1997; Akeson et al. 1996). They are composed of an aluminum mirror and a grid of gold- 
plated tungsten wires of 25 /zm diameter and spaced at an interval of 125 /im. The inside diameter 
of the grid is roughly 16 cm, some 25 times bigger than the incident Gaussian beam at 1.3 mm 
(Akeson 1997). 

In this section, we will compare data obtained in the calibration of these polarizers at 1.3 mm 
with the model calculated earlier. In the experimental set-up, the incident beam is composed of 
radiation emanating from a hot load (absorber at room temperature) polarized along the vertical 
axis and a cold load (absorber in liquid nitrogen) along the horizontal axis. The beam is incident 
on the polarizers at an angle of Xi = 34 deg. with the grid rotated by <p g = ±50.3 deg. relative to 
the vertical, all in the coordinate system of the laboratory (coordinates (it,v,iu) of Figure 1). These 
values can be inserted in the appropriate equations of our earlier analysis and used to test our model 
against the experimental data. 

The calibration consists in using our model to map out the actual distance between the grid 
and the mirror as the latter is moved with a micro-positioner which is part of the assembly. When 
this is done, this distance can then be precisely adjusted to ^ in order to use the polarizer as 
a reflecting quarter- wave plate for polarimetry measurements. The grid rotation angle tp g must 
also be calibrated so that it can be set to the proper value that will allow the transformation 
of incident linear polarization to circular polarization. (This condition is met for ft — ±ay (or 
tan(^) = co ^ y ), as can be asserted from our earlier discussion of the principal axes of a grid in 
section 3.3; this gives <p g = ±50.3 deg as quoted above). 

Figure 4 shows the results obtained from such measurements (of the reflected polarized intensity 
along the horizontal -tZ-axis in the laboratory coordinate system) made on antenna #6 of the OVRO 
array at a frequency of 232.037 GHz when the separation between the grid and the mirror is varied 
through a range of a several hundreds of microns. Accompanying the data points is a least square fit 
of our model (solid curve) with no free parameters as far as the grid is concerned, only the hot and 
cold load levels and the offset in the mirror-grid separation were allowed to be fitted. The agreement 
is very good. The main shortcoming of the fit is at a backshort position of roughly 900 /xm where 
a resonance is evident from the data. The model also shows a resonance at the same position but 
the fit is not perfect. This feature is caused by the small amount of unwanted transmission from 
the component of the incident electric field aligned with one of the principal axes (pi) which gets 
trapped between the grid and the mirror. 

Before we try to explain the differences in width and shape of the resonance, we would first like 
to show two ways by which it can be suppressed (such a response from the polarizer is a nuisance 
when trying to calibrate it and should be avoided). 
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First, reducing the amount of unwanted transmission through the grid would certainly have a 
damping effect on the resonance. We have shown how to do just that in section 3.4 when defining a 
set of optimum parameters for a grid, so using equations (64) and (65) we find for our application 
a ~ 24 /im and d ~ 148 /mi. Figure 5 shows a comparison of the simulated responses for the 
polarizer tested at OVRO and our optimized polarizer. As can be seen, any sign of the resonance 
has disappeared in the latter. 

Another way of avoiding the resonance, while still using the same original grid with a = 12.5 /im 
and d = 125 /im, is to replace the mirror with another grid (Young 1997) and rotate both of them in 
such a way that the projected orientation of their wires in the plane of the incident field are aligned 
with a principal axis. By this we mean that the two grids have their angle of rotation specified 
by tan(y> 9 ) = ( or *9 = 503de g) and tan (^) = (or ^ = -50.3 deg) respectively. 

This would ensure that the unwanted transmitted field from the first grid would almost be entirely 
transmitted through the second grid, therefore getting rid of the resonance. 

Obviously, trying to solve for such a configuration using Maxwell’s equations would he a 
formidable task. We can however use our scattering matrix model developed in section 3.3. We 
then have to define two matrices, one for each grid, and solve the problem for cases where they are 
separated by a given distance while terminating the last grid by the line characteristic impedance 
Z 0 {E£ = Ef = 0). We simulated the response predicted for such an arrangement of grids and got 
results that are practically identical to those presented in Figure 5 for the optimum polarizer. 


5. Grids and beams of radiation. 

Until now we have restricted our analysis to cases where the dimensions of the grid (or the 
assembly) and the extent of the incident wave were assumed to be infinite. These simplifications 
were necessary in order to allow us to have a chance at a solution, as the reality of finite sizes brings 
severe difficulties in the analysis. It would however seem reasonable to suppose that if the incoming 
excitation can be properly represented by a beam of radiation which is of a size a few times smaller 
that the actual dimensions of the assembly that the results obtained with our analysis should still 
be valid. Indeed, one could argue that the incident beam should only induce currents in the vicinity 
of the area where it impinges on the assembly. There should therefore be little to no differences in 
its response whether it is infinite or not. 

Although we believe this argument to be a reasonable one, we will show that the characteristics 
of the incoming radiation can be important in some cases. We will in fact argue that it can explain 
the discrepancies in the width and shape of the resonance observed in the response of the reflecting 
polarizer presented in section 4.3 (see Figure 4). 

As a starting point, let’s take note that we can always mathematically express a beam of 
radiation i? 0 (r) as a summation of plane waves with different amplitude (and phase) and k vectors. 
For example, in the laboratory system of coordinates (u,u,u/), the electric field along the u-axis 
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Eii r ,t) — a " E 0 {t) exp (jui n t.) can be expressed by (using its Fourier transform in (w, k)-space and 
assuming the beam to be monochromatic at lj 0 = kc): 


E 0 ( k,w) 


a" f°° 

,r-T 4 / du)d 3 kE 0 (k,w)exp[-j {k-r -wt)] 
(27tJ J - oo 


/ oo 

d 3 r / £ , 0 (r / ) exp (jk * r') 

-oo 


where S ( x ) is Dirac’s delta distribution. 


(81) 

(82) 


Let’s now assume that the incident radiation can be satisfactorily modeled using a circular 
Gaussian beam with a beam waist W Q and a Rayleigh range z R = We also know that the 

resonance will occur for a grid-mirror separation of h ~ jf- for each spectral component, where k z is 
the projection of the wave vector along the z-axis perpendicular to the reflecting polarizer (in the 
coordinates system of the grid of Figure 1 ). From this we can express the width of the resonance 
Ah as a function of k z and A k z the spectrum extent along the same axis: 


Ah ~ 7 r 


A k z 

"fcP 


We need to find an expression for A k z and this can be done as follows. Using the wave 
uncertainty relation, we can evaluate the spectrum extent in the laboratory system of axes as: 


Afc u 

~ W~ 1 

— rr O 

A.k v 

~ w~ l 

— rr 0 

-V 

<1 

~ z~ l 
- Z R . 


Transforming these in the coordinate system of the grid and inserting the result in the expression 
for the width of the resonance we find: 


When we use the corresponding values of the different parameters appearing in equation (83) 
for the case of the reflecting polarizer discussed in the last section (A = 1.3 mm, 7 = 0.83 and 
W 0 — 3 mm) we find Ah ~ 37 fim. If we take into account the finite bandwidth of the OVRO 
receivers (1 GHz), a similar exercise shows that at most only a few microns need to be added to 
the previous estimate. Although these numbers represent only a rough calculation of what could 
be expected, they nevertheless tell us that there will be a significant broadening of the resonant 
feature. 
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We will not try to produce a perfect fit to the data obtained at OVRO for this would require 
extensive modeling of our experimental set-up and therefore bring us to a level of complexity that 
we do not wish to tackle at this time. But using equations (81) and (82) and applying the result of 
our analysis of the reflecting polarizer for every spectral component hence calculated, we can get a 
better idea of the phenomenon considered here. 

We have done this and the result is shown in Figure 6 where we present the result of a simulation 
of the effect of a Gaussian beam on the width and shape of the resonance exhibited by a reflecting 
polarizer of the kind discussed in section 4.3. The beam is converging with its waist situated some 
10 cm “behind” the polarizer and the integrated power (over a beam width) is measured some 
distance away from the assembly in the far-field. Although there still remain some differences, this 
simulation shares a lot of the same features observed experimentally. 

We believe that simulations like this one along with our earlier calculations provide convincing 
and compelling evidences for the importance of appropriately taking into account the nature of the 
incident radiation in the analysis of similar systems. 


6. Conclusion. 

In this paper, a general solution for the analysis of polarizing grids was presented; it is valid for 
arbitrary angles of incidence and of grid rotation. With it and the scattering matrix representation 
that derives from it, basically any configuration or system of grids can be analyzed as long as some 
assumptions concerning the wire radius and spacing are respected (A > 40a and d > 4a). This is 
not a severe restriction as most grid currently available satisfy those conditions, we refer the reader 
to Chambers et al. (1986, 1988) for cases where a larger size of wire is needed. Our analysis also 
allowed us to define a set of optimum values for both the wire radius and spacing as specified by 
the following equations: 


L(1 -a 2 ) 4 7r 7 <7Zj 

d ~ 2tt a . 

We provided an analysis of the effects that two types of random errors can have on the per- 
formance of a grid. It was shown that errors in the wire spacing were the most important and 
could have some impact on the amount of unwanted polarization transmitted through a grid. In 
that respect, our model showed to be in good agreement with the experimental results of Shapiro 
& Bloemhof (1990). 

Comparisons with experimental data obtained in the calibration of a reflecting polarizer used 
at the OVRO were also presented and predictions from our model are in good agreement with it. 


The only discrepancies appeared in the nature of a resonance, more precisely its width. But we 
have shown that it could be accounted for by including in the analysis a proper treatment of the 
effects of the nature of the incident radiation on the response of the polarizer. 

We are grateful to J. B. Shapiro and E. E. Bloemhof for their permission to use their previously 
published experimental results. We wish to thank the staff of the Owens Valley Radio Observatory, 
and O. P. Lay for numerous discussions and suggestions. The Owens Valley Radio Observatory is 
funded by the National Science Foundation under Contract No. AST 96-13717 and the polarimetry 
project at OVRO through NASA grant NAG5-4462. M. H. work was supported in part by a grants 
from FCAR and the departement de physique de FUniversite de Montreal. 
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A. Generalization of the principal axes. 

Referring to equations (49) for the scattering matrix, one first realizes that such a symmetry 
in its components will always be seen when representing an arbitrary polarizing grid. Only the 
functions which define R uu y T uu , .. . will change. Further, when trying to reduce the scattering 
matrix to a form similar to equation (56) it is only necessary to concentrate on only one of the 
two blocks (each appearing twice) present in equation (49). For example, if one diagonalizes the 
block composed of the reflection coefficients then the transmission block is also diagonalized and 
vice-versa. In obtaining the results which follow, we have worked with the reflection block appearing 
in the upper left and lower right of equation (49). We will now show that the orientation of the 
principal axes is determined by the symmetry of the grid and can be deduced using the formalism 
of group theory. 

As seen by the incident wave, the grid has a symmetry which can be expressed by a represen- 
tation of the point group C<i v The four covering involved are: the identity (E). a rotation by 7r 
about the w-axis (C 2 ), a reflection (<r v ) across a plane defined by the re-axis and and an axis defined 
by the projection of the direction of the wires in the plane of the incident field and finally another 
reflection (< 7 f v ) across a plane perpendicular to the previous one (and to the plane of the incident 
field). Upon studying the character table of this group (see Tinkham (1964, p. 325)) and the effect 
of the above operations on the two possible states of linear polarization (along the u and t>-axes) 
we find that only two non degenerate irreducible representations ( B\ and B 2 in Tinkham (1964, p. 
325)) will be realized. For each of these there will exist one eigenvector, each corresponding to a 
given principal axis. These can be deduced by constructing the appropriate symmetry coordinates 
(Wilson et al. 1955) which turn out to be the two principal axes pi and P 2 previously defined in 
section 3.3.1. 

Since this result was obtained with the use of group theory, it is perfectly general and indepen- 
dent of any approximations that can be used in dealing with a polarizing grid. 
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B. List of symbols. 

a wire radius, 

A s scattering vector potential, 
c speed of light in free space, 
d wire spacing, 

Ej,E 5> E^ incident, scattered and total electric field, 

F constant (= ^ a;Q ) , 

F s scattering vector potential, 
h grid-mirror separation (reflecting polarizer), 

H s scattered magnetic field, 

Hn Hankel function of the second kind of order n, 

J current density vector, 

J n Bessel function of the first kind of order n, 
k wave vector of the incident wave (|k| = k = %L), 
k' = kVT^a*, 

K total surface current density vector, 

K X ,K & mean longitudinal and azimuthal surface current densities, 

Pi ^ p 2 principal axes of a grid, 

K x ,R y ,R z reflection coefficients in the system of coordinates of the grid (see Figure 1), 

R u , R v reflection coefficients in the system of coordinates of the laboratory (see Figures 1 
and 2), 

R\\,R± reflection coefficients along the principal axes of a grid, 

Rte , Rtm transverse electric and transverse magnetic reflection coefficients of the mirror, 

S scattering matrix, 

T x ,T y ,T z transmission coefficients in the system of coordinates of the grid (see Figure 1), 

T U ,T V transmission coefficients in the system of coordinates of the laboratory (see Figures 1 
and 2), 
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7j|,7j transmission coefficients along the principal axes of a grid, 
W 0 <ZR beam waist and Rayleigh range of a circular Gaussian beam, 


Z impedance matrix, 

Z a , impedance of free space (= 

Z p , Z n grid impedance along the principal axes (as defined with the impedance matrix), 
Z s surface impedance of the wires (= (1 + j) \J^)i 

Z\\, Z L grid impedance along the principal axes (as defined with the scattering matrix), 

a projection of the normalized wave vector on the x-axis 
(= sin (xi) sin (<£«,)), 


j3 projection of the normalized wave vector on the y-axis 
(= sin (xi) cos (<^ 9 )), 

•y projection of the normalized wave vector on the z-axis (= cos (x»))i 
a',/?', 7' projection of the normalized incident field on the x, y and 2-axes, 
a" ,(3" projection of the normalized incident field on the u and u-axes, 


6 (x) Dirac’s delta distribution, 

E 0 , fx a permittivity and permeability of free space, 
t 7, £ random errors in wire radius and spacing, 

A wavelength, 
a wire conductivity, 

(p Xi angle of grid rotation and angle of incidence, 
Euler’s constant (^ 0.577215), 
u angular frequency of radiation. 
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Fig. 1. — Coordinates system for the study of a polarizing grid or a single wire. The wave vector 
k of the incident radiation is aligned with the re-axis, the u, x and y-axes are in the plane of the 
page, the w and z-axes are in the plane perpendicular to the u-axis (into the page) and the wires 
are parallel to the zy-plane. We refer to the (u,u, w) and (x,y, z) systems as the laboratory and 
grid coordinates respectively. 
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Fig. 2. — Definition of the system of coordinates {u\ v\ w f ) for the reflected wave in relation to the 
(u, v,w) system of the incident /transmitted waves introduced earlier in Figure 1. The u ! and u-axes 
are one and the same and are pointing out of the page. The direction of propagation of the reflected 
wave is along the negative w'-axis. 
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Fig. 3.— Curves of predicted values for the cross- polarized transmittance plotted against experi- 
mental data from Shapiro & Bloemhof (1990) and Bloemhof (1998). The three grids have a random 
error (1-a) in wire positioning of 5%, 16% and 37% with mean distance between wires of 103pm, 
109pm and 114 pm respectively; they all have a wire radius of 12.5 pm. 
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Pig 5 — Comparison of the predicted results obtained for the polarizer tested at OVRO (solid 
curve) and our optimized polarizer (broken curve). The resonance is not present on the optimized 
polarizer’s response. The two curves are plotted with a small vertical offset between them. The 
intensity is in arbitrary units. 








